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Abstract 

We consider a (2 + l)-dimensional SU(A) lattice gauge theory in an axial 
gauge with the link field U\ set equal to one. The term in the Hamilto- 
nian containing the square of the electric field in the 1-direction is non-local. 
Despite this non-locality, we show that weak-coupling perturbation theory 
in this term gives a finite vacuum-energy density to second order, and sug- 
gest that this property holds to all orders. Heavy quarks are confined, the 
spectrum is gapped, and the space-like Wilson loop has area decay. 



1 Introduction 



The central problem of QCD is confinement. It is not enough to prove that 
lattice gauge theories have a confining phase - which is evident from strong- 
coupling expansions. It is necessary to see that the color is confined at 
arbitrarily weak lattice coupling. 

We find here that (2 + l)-dimensional lattice gauge theories confine for 
any dimensionless bare coupling. The technique used is a weak-coupling 
expansion in an anisotropic lattice gauge theory in an axial gauge. This is not 
the standard expansion utilizing Feynman diagrams. Though the coupling 
constant of a (2 + l)-dimensional gauge theory is not infinitely renormalized 
in (2 + 1) dimensions, the dimensionless bare coupling on the lattice must 
vanish in the continuum limit. This is why a weak-coupling analysis is useful, 
even for this case. The dependence we find of the string tension and the mass 
gap on the coupling constant does not agree with conventional wisdom - our 
results for these physical quantities don't behave as anticipated, as the lattice 
spacing is taken to zero - but they are not zero. 

The first analytic demonstration of confinement of heavy sources in (2 + 
l)-dimensional gauge theories was given by Polyakov for lattice compact 
QED pp, and later for the Georgi-Glashow model [2]. The latter model is 
interesting in that color charges disappear completely from the spectrum. 
This is, however, different from the sort of confinement we expect for QCD, 
in that matter fields play an important role. Feynman argued that (2 + 1)- 
dimensional QCD is confining 3J. Unfortunately, the orbit-space distance 
estimates in Feynman's paper are incorrect. Nonetheless, his basic claim, 
that the diameter of gauge-orbit space of (2 + l)-dimensional SU(iV) Yang- 
Mills theory for small magnetic energy is finite, appears to be correct [3]. New 
nonperturbative methods which do not require a lattice have been derived 
by Karabali, Kim and Nair [3] (one of their formulations of the Hamiltonian 
has been obtained a different way in reference (HI). 

There is one special set of assumptions we use to derive our results; the 
(l + l)-dimensional non-Abelian nonlinear sigma models (without topological 
terms) have a mass gap, exponentially decaying correlation functions, and 
their vacuum expectation values of local operators exhibit clustering [7j, [S], 
jH]. Though no rigorous proof of these properties exists, we think that the 
evidence in their favor is overwhelming. 

Our basic strategy is to write the lattice version of the Hamiltonian as 
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the sum of two terms, namely 



H Q = I d 2 x { ^Tt£ 2 2 + ^TtB 2 



and 



e 2 e 2 



—V=— d z xTr£( 
2 2' 1 



where £j are the components of the electric field conjugate to the gauge 
field Aj, [£j(x), A k (y)} = i5jkS 2 (x — y) and B = \[d\ — L4i,d 2 — iA 2 ] is 
the single space component of the magnetic field. We then pick the gauge 
A\ = 0. When this is done on the lattice, H is a set of decoupled chi- 
ral SU(iV) xSU(iV) nonlinear sigma models for which the S'-matrix and the 
spectrum are known. The quantity V is non-local, but we show that per- 
turbation theory in this term is sensible to second order. The vacuum state 
in this perturbation series confines fundamental color charges. Our splitting 
of the Hamiltonian is not explicitly rotationally invariant, but if the method 
works to all orders of perturbation theory, rotational invariance should be 
restored. 

Let us review axial gauges in the continuum. If the SU(iV)-Lie-algebra- 
valued gauge field A\ is set to zero, then Gauss's law may be integrated to 
obtain 



x i d-i 

£ x (x) = - dy x Yyd j -iA j {y\x 2 ,...x d - 1 ) 1 £ j {y\x 2 1 ... 1 x d -" 

3=2 



dy 1 D ± (y 1 ,x L )-£ ± (yW) , (1.1) 

where the dimension of space is d — 1, x 1 - = (x 2 , . . . , x^" 1 ), and D± are the 
last d — 2 components of the covariant derivative in the adjoint representation 
(Z>2, • • • , Dd-i). The term in the Hamiltonian 

Y V = J d^x^T^l (1.2) 

must have a vacuum expectation value proportional to the volume of (d— 1)- 
dimensional space, if the theory is to be sensible. As discussed by Mandel- 
stam PU] this means that the quantity 

K(y\ z\ x ± ) = <0|Tr D ± (y\ x^ -£ ± (y l , x L ) D ± (z\ x^ -^(z 1 , x x )\0) 
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must have the property that J dy 1 dz l K(y 1 , z 1 , x x ) does not diverge with the 
spatial volume. One might think that if K(y x , z 1 , x ) falls off sufficiently fast 
with ly 1 — z 1 \, the problem can be ameliorated. Rapid fall-off of K, however, 
is not enough. Even if the fall-off is exponential, the result may diverge as 
(L 1 ) 2 where L l is the range of x 1 . Mandelstam recognized that the residual 
gauge invariance, remaining after solving for E\ in 1)1.1)1 . namely 

J dx l D ± ■ £ ± * = , (1.3) 

must also be satisfied by the vacuum. Without both ([1.3)1 and and sufficiently 
rapid decay of K(x 1 ,y l ,x ± ), any conjecture for the vacuum may have an 
unacceptable infrared-divergent energy, coming from (|1.2p . Fortunately, we 
find that in our perturbation scheme, both the unperturbed vacuum energy 
and the first two corrections in our weak-coupling expansion obey the lattice 
versions of both the rapid-decay criterion and ([1.3)1 . 

2 The lattice gauge Hamiltonian 

The purpose of this section is to establish our definitions and conventions. 
It is not an introduction to the Hamiltonian SU(iV) gauge theory. Such 
introductions can be found in the review article by Kogut and in the book 
by Creutz [TT] . 

Consider a lattice of sites x of size L l x L 2 , with sites x whose coordinates 
are x l and x 2 . We require that x 1 /a and x 2 /a are integers, where a is the 
lattice spacing. There are 2 space directions, labeled j — 1, 2. Each link is a 
pair x, j, and joins the site x to x + ja, where j is a unit vector in the j th 
direction. 

We introduce basis vectors or generators t a , a = 1, . . . , N 2 — 1, of the 
Lie algebra of SU(iV). Sometimes we use Roman letters for the index, e.g. 
we may write t\, rather than t a (the purpose of using different alphabets 
is to distinguish between coordinate indices on the SU(iV) manifold and 
tangent-space vectors). The generators are defined to be orthonormal, so that 
Tr t a t/3 = 5 a p. The structure coefficients of the Lie algebra, /2g, a, (3, 7 = 
1, . . . , iV 2 — 1, are, as usual, the complex numbers defined by [t a , tp] = if^pty 
The identity matrix will be denoted by 11. 

The Hamiltonian lattice gauge theory is usually formulated in temporal 
gauge A = 0. The basic degrees of freedom, before any further gauge fixing, 
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are elements of the group SU(iV) in the fundamental (N x N) -dimensional 
matrix representation Uj(x) G SU(iV) at each link x, j. In addition, there are 
the N 2 — 1 electric-field operators at each link lj(x) b , b = 1, . . . , N 2 — 1. The 
electric-field operators are self-adjoint by construction. The commutation 
relations on the lattice are 

[lj(x) b ,l k (y)c] = i$xy$jk fbc h{ x )d , 

[lj(x) b ,U k (y)] = -S xy 5 jk t b Uj(x) , (2.1) 

all others zero. In the Schrodinger representation, with the components of 
Uj(x) taken to be c- numbers, the latter of (|2.1j) becomes 

lj{x)bU k {y) = SxySjk h Uj(x) , 

The lattice Hamiltonian is 

2 iV 2 -l 2 i 

^ = EEE |[^] 2 -Ei^ [TrU 12 (x)+TrU 21 (x)} , (2.2) 

x j=l 6=1 x ^° 

where 

Uj k (x) = Uj(x)U k (x + ja)Uj(x + ka)^U k (xy , 

and the bare coupling constant go is dimensionless. Note that the coefficient 
of the kinetic term can be written in terms of the continuum coupling con- 
stant e, namely ^q/(2o) = e 2 /2. It is for this reason that hadron masses and 
the string tension evaluated in lattice strong-coupling expansions all scale 
sensibly with e, in (2 + 1) dimensions. 

We denote the adjoint representation of the SU(iV) gauge field by 71: 

N 2 -l 

K b % = Ut b U ] . 

c=l 

The matrix 1Z lies in the group SU(iV)/Z7v. This is a special orthogonal 
matrix n T TZ = 1, det K = 1, and ^\]{N)/Z N is a subgroup of SO(iV 2 - 1). 

Schrodinger wave functions are complex-valued functions of all the link 
degrees of freedom Uj(x). Physical wave functions \I/ ({£/}) satisfy Gauss' law 

2 

(V ■ l)(x) b *({£/}) = Pi h(*)] b *({U}) = , (2.3) 
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where 



N 2 -l 

[Vjljix)^ = lj(x) b - n i( x ' C hi? - ka) c . (2.4) 

c=l 

Sometimes it is useful to introduce color charge operators at lattice sites, 
denoted by q(x) b , which satisfy 

[q(x) b , q(y) c ] = ifb c 6 xy q(x) a . (2.5) 

In the presence of charges, Gauss's law becomes 

[(V.l)(x) b -q(x) b \ *({[/}) = 0. (2-6) 

Henceforth, we will drop the explicit summation symbol for repeated 
group indices and adopt the Einstein summation convention. Sometimes we 
omit the lattice site or link labels, provided no confusion should be caused 
by such omissions. 

There is a natural geometric interpretation of the electric-field operator. 
The Maurer-Cartan vector e a a , on the manifold of SU(iV) defined by 

e a a t a = -iU~ l d a U , 

is given explicitly by 

in canonical coordinates A a , a — 1, . . . iV 2 — 1, defined by U = e~ tA ' t , and 
d a = d/dA a . The coordinates A are related to the continuum gauge field A 
by A = a A. The matrix e is nonsingular (including at A a = 0). One may 
view e a a as the linear map from the group manifold to the tangent space; 
this is a particular choice of the vielbein, and in this case there is torsion. 
The electric-field operators are given by 

l a = -i{e-\ a d a . 
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3 The axial gauge on a cylinder 



By fixing an axial gauge, we will find the gauge- invariant degrees of freedom, 
up to coordinate singularities of measure zero. Such gauge fixings have been 
discussed many years ago, both in the continuum ITSj and on a lattice [13], in 
the path-integral approach to gauge theories. The advantage of working with 
the Hamiltonian instead of the path integral is that unphysical components 
of the gauge fields may be more easily eliminated using Gauss's law £[4J (this 
could also be done in a transfer matrix formalism). 

We choose space to be a lattice cylinder of size L 1 x L 2 , with periodic 
boundary conditions in the 2-direction only. This means that for any function 
fix 1 , x 2 ) of lattice sites fix 1 , x 2 + L 2 ) = fix 1 , x 2 ). We take components of x 
to have the values x 1 = 0, a, 2a, . . . , L 1 , and x 2 = 0, a, 2a, . . . , L 2 — a. Gauss's 
law is still given by ()2.3|) . provided (|2.4j) is modified to 

£>i/i(x) = (1 - 5 x i L i)hix) - (1 - 5 x i o^iix 1 - a,x 2 )k{x l - a,x 2 ) , 
T>%li{x) = h{x) — Hzix i x 2 — a^ix 1 , x 2 — a) , (3.1) 

to take into account points on the boundary. 

We gauge-fix the links in the 1-direction by ^(x 1 ,^ 2 ) = 11 everywhere 
and use ()2.3|) and ()3.1|) to write 

l 1 {x\x 2 ) = -Y,(P 2 -h){y\x 2 ). (3.2) 

There is some non-Abelian gauge invariance remaining, namely that 

L 1 

T(x 2 )^ = ^(D 2 ■ / 2 )(x 1 ,a; 2 )^ = . (3.3) 

We split the Hamiltonian into two terms H = H + kV, where eventually 
we set k = |r 

- ^2— r [TiU 2 {x\x 2 ^U2ix 1 + a,x 2 ) + c.c.}\, (3.4) 

xi=0 29 ° a ) 
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and 



L 2 -a L 1 



" = EE 

x 2 =0 x!=0 



2,1=0 



(3.5) 



It will be important for the discussion in the next section that the con- 
straint (|3.3j) allows us to replace (|3.5|) by 



x 2 =0 x!=0 



n T 



"£(V 2 -l 2 )(y\ 

j/l=0 



.z- L =a; 1 +a 



. (3.6) 



We have assumed until now that no charges are present. If a quark is 
placed at site u, then ()2.6|) may be solved to give 



hix 1 ^ 2 ) = q(u\u 2 ) 5 x i> u i 5 X 2 U 2 - £(£> 2 - l 2 )(y l , 

3,1=0 

The remaining gauge invariance is 



x 



(3.7) 



£(£>2 ■ Z 2 )(a a ,2 2 ) - q(u l ,u 2 )5 x 2 u 2 



.xi=0 



* = 



(3.8) 



4 Confinement at leading order 

The splitting (J3.4)) . (J3.5)) is not 90° rotation invariant. Nonetheless, if pertur- 
bation theory in V makes sense, this rotation invariance should be restored 
at sufficiently high orders. Notice that H is a set of decoupled (1 + 1)- 
dimensional lattice chiral non-linear sigma models, with global symmetry 
SU(iV)i x SU(iV)fl, plus an extra term at the boundary x 1 = (this is 
a sum of unitary matrix- model Hamiltonians) . The on-shell properties of 
these sigma models have been completely determined; the Bethe Ansatz jS] 
and analytic 5*-matrix theory determine the spectrum in the renormalized 
continuum limit. 

Let us briefly describe the particles of the chiral SU(iV) l x SU(A^)^ model. 
There are fundamental particles with mass m x transforming as the fully an- 
tisymmetric tensor representation of SU(N)l x SU(iV)#. The particles are 
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labeled by a quantum number n = 1, . . . N — 1. The particle with n > 1 is 
a bound state of n fundamental particles. We may regard the bound state 
of n particles as a bound state of N — n antiparticles. There is no singlet in 
the one-particle spectrum (which would correspond to n = N) . The particles 
have masses 



m. 



mi 



sin^| 



n 



N 



sin 



N 



where the mass gap nil is of the form 
a \ J a 



mi 



/ 2 \^ Ki 

(e a) 2 exp — — + 



(4.1) 



where C is a non- universal constant, K\ = —\ and K2 = 47r are determined 
from the one- and two-loop coefficients of the chiral-model beta function [TH] . 
respectively, and the corrections are non- universal. 

Suppose that there are no charges present. The remaining gauge invari- 
ance ()3.3j) means that we impose on the states ^ of the chiral model at x 2 
the constraints 



Riix 1 , x 2 - a)l 2 (x l , x 2 - a)* = V] hfa 1 , x 2 )^ 



(4.2) 



The meaning of (|4.2j) is that if the state of the chiral model at some particular 
x 2 transforms as a vector with some set of weights under SXJ(N)l, then the 
state of the chiral model at x 2 + a transforms the same way under SU (N)r. 

For the ground state which is a product of chiral model ground 

states, each side of (|4.2|) is automatically zero; for the Hohenberg-Mermin- 
Wagner theorem guarantees that it is a singlet under both the left global 
SU (N)l and under the right global SU(iV)# invariances. 

The leading-order Hamiltonian describes a theory which confines funda- 
mental charges separated in the 2-direction. We will show that this is true by 
two different lines of reasoning. The first proof is more in line with the way 
people usually think about phenomena in gauge theories. The second proof 
is a direct utilization of the concepts we have used in the previous section 
and this one. 

Here is the first proof: suppose that a quark is placed at w 1 ,^ 2 and an 



anti-quark at m 1 ,?/ 2 ^> u 2 . Gauge-invariant states are of the form 
\C >= A(u\ y y Yl UQink) B(u\u 2 y |0 > , 



(4.3) 



linkeC 
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for some path C of links joing the quark to the anti-quark, whose creation op- 
erators are A* and B\ respectively. The lowest-energy state in the presence of 
the sources is a superposition of such states. The Hohenberg-Mermin- Wagner 
theorem states that for a Hamiltonian with a global continuous symmetry, 
there is no spontaneous symmetry breaking. In the unperturbed vacuum, 
therefore, < 0[ C/21 >= 0. This means that the action of U 2 on the chiral- 
sigma-model ground state produces a superposition of excited states only. 
Thus the expectation value of Hq in any state (|Oj) must be bounded below 
by the gap times the separation of the fundamental charges, i.e. 

< C\H \C> > - u 2 \ , (4.4) 

a 

which means that there is confinement of fundamental charges, with string 
tension m\j 'a. We call this phenomenon "vertical confinement", because 
confinement occurs in the 2-direction. 

Now for the second proof: the constraint (|3.8jl has the form 

L 1 

R 2 (x 1 , x 2 — a)/ 2 (x 1 , x 2 — a)ty — q^u 1 , u 2 )S u 2 x 2^/ + q{u x , u 2 )5 u 2 y 2^/ 

x x =0 

L 1 

= ^Z 2 (x 1 ,x 2 )*. (4.5) 

This tells us that if the chiral model at u 2 — a is in an SU(iV) ^ singlet state 
(such as the vacuum), then the chiral model at u 2 cannot be in an SU(A^)^ 
singlet. Thus the chiral model at u 2 is in an excited state. By continuing to 
use (J4.5|) we conclude that all the chiral models for x 2 satisfying u 2 < x 2 < y 2 
are excited. In this way, we obtain the same result for the vertical string 
tension as that given above. 

A rectangular Wilson loop of size S\ x 5* 2 is 

A(S 1 x S 2 ) = Tr W(x\ x 2 ; S 2 ) ] W(x l + S u x 2 ; S 2 ) , (4.6) 

in our gauge, where 

W{x\ x 2 - S 2 ) = U 2 (x\ x 2 + S 2 )--- U 2 (x\ x 1 ) . 

Correlation functions of U 2 decay exponentially. We expect that for large Si 

<0\[U 2 (x\x 2 y} a b U 2 (x 1 + S lj x 2 ) c d \0>~ D b a iex V (-m 1 S 1 ) . 
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The Wilson loop expectation value is a product of S 2 /a such correlation 
functions, and therefore 

< 0\A(S 1 x S 2 )\0 > ^ exp (-^S^ • (4.7) 

This is an area law, with the same string tension mi/ a found above. 

There is no "horizontal confinement" - that is, there is no confinement in 
the 1-direction - yet. Horizontal confinement will only appear if the pertur- 
bation kV = <7oVy(2a) is taken into account. This is because the constraint 
consistent with the presence of a quark at x l ,x 2 and an anti-quark at y l ,x 2 
with y 1 ^> x 1 is ()3.3|) . which is satisfied by the unperturbed vacuum. Thus, 
if V is neglected, there is no force between a quark- ant i- quark in the 1- 
direction. The appearance of horizontal confinement in perturbation theory 
will be demonstrated in Section 6. 



5 Weak-coupling perturbation theory and in- 
frared finiteness 

If L 1 and L 2 are kept finite, the spectrum of the Hamiltonian H + kV is 
purely discrete. Let us consider this spectrum to second order in Rayleigh- 
Schrodinger perturbation theory: 

E n = E^ + KE^ + . 2 E^ + ... , 

where 

I <r H/ (0) lT/|vR (0) > |2 

4" =< *Sf |v|»Sf >, *p = - £ l< y^ ) >l . ■ ■ ■ , (5.D 



and \^n^ > are the eigenvectors of H with eigenvalues En\ The purpose 
of this section is to show that the corrections to the vacuum energy E^ 
and E^ are proportional to L 1 x L 2 . This happens for two reasons: 1) the 
ground state of the chiral model is "disordered" , i.e. two-point functions fall 
off exponentially, and 2) the unperturbed vacuum is a singlet, simplifies the 
form of V acting on this vacuum to (|3.6|) . Our philosophy is close to that of 
Mandelstam [TQJ in this regard. 
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The first correction to the vacuum energy is 

i#> = - < E E PWvV)F 

s l=0 y 1 =0 z 1 =x 1 +a 

Correlation functions of Z 2 and -R 2 ^ 2 must decay exponentially with the dis- 
tance, and therefore this quantity will have the form 

fl«~-L a EE E e-^^l, (5.2) 

The dominant contribution to this expression comes from y 1 « z . Since 
y 1 < a; 1 < z 1 , £7q is proportional to the volume L 1 L 2 . 

Next we sketch the proof that the second-order correction to the vacuum 
energy also scales linearly with the volume. Notice that the coefficient of 
each energy denominator in the second correction (|5.1j) is non-positive. Thus 



I A) 



< ^Ei<*fW>i ! 

1 m^O 

= -1 [<vl/( )|^| v]> (0) > _ (< ^(0)^^(0) >)2 j ^ (5 3) 

The connected vacuum expectation value on the right-hand-side of ()5.3|) has 
the following form: 

L 2 L ' 



S 2) i < ^EE^-v 2 ), 



where 



cv,»V) = EE E E E [<» 



(0)| 

I 



r y 1 =0 z 1 =x 1 +a u 1 =0 v 1 =w 1 +a 
.1 ~2\7Vi ; / 1 2\ 



x V 2 l 2 (y L ,x 2 yV 2 l 2 (z\x 2 
x V 2 l 2 (u\x 2 + ia) T V 2 l 2 ( y v\x 2 + ra)|^ 0) > 
- <^ 0) |I? 2 / 2 ( 2/ 1 ,x 2 ) T D 2 / 2 (z 1 ,x 2 )|^> 

x< ^ 0) |I? 2 / 2 ( M 1 ,x 2 + ra) T I? 2 / 2 (t; 1 ,x 2 + ra)|^[ ) 0) >1 , (5.4) 
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and where r = 0,±1. The chiral model is a massive local quantum field 
theory, so that vacuum correlation functions must cluster for the dominant 
part of the summations in (|5.4|) . Therefore this expression is approximated 
well by 



c(x\wV) « J2E E E E [<^> 2 ^v 2 f 

r y 1 =0 z 1 =x 1 +a u 1 =0 v 1 =w 1 +a 

x ^(liS^ + ra)!^^ 

x <^\V 2 l 2 (z\x 2 ) T V 2 l 2 (v\x 2 + m)\^ ] > 

x <^f\V 2 l 2 {z\x 2 ) T V 2 l 2 {u\x 2 + ia)\^ ] >\ . (5.5) 

By using ()3.3|) we can write each term of ()5.5|) as something which vanishes 
exponentially away from x 1 = w 2 . For example, consider the first factor of 
the first term: 

X 1 z 1 

First Factor = E < ^o ) \ D 2h{y\x 2 ) T V 2 l 2 {u\x 2 + ra)|^ 0) > . (5.6) 
If x 1 < z 1 , we may write this as 

x 1 L 1 

F.F. = -J2 E <^t\v 2 l 2 {y\x 2 ) T V 2 l 2 {u\x 2 + ia)\yf>, (5.7) 

y 1 =0 ■u 1 =z 1 +a 

and we see that this expression is finite as L 1 —>■ oo. On the other hand, if 



x 



1 > 2 1 , we rewrite (15. 6 J) as 



F.F. = - ^ ]T <^f\V 2 l 2 {y\x 2 ) T V 2 l 2 {u\x 2 + ia)\^ ) >, (5.8) 

j/ 1 =x 1 +a n 1 =0 

and reach the same conclusion. Since each factor of each term behaves this 
way, we can conclude that the second-order correction to the vacuum energy 
can increase at most linearly with L 1 . 

Infrared finiteness of the vacuum energy to first and second order in per- 
turbation theory inspires confidence that it should hold to all orders. The 
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main complication beyond the second order is the lack of non-positivity or 
non-negativity of products of matrix elements. We believe that careful ap- 
plication of the linked-cluster expansion, assuming clustering in the chiral 
sigma model, can provide a proof to all orders. 

6 Horizontal confinement 

In Section 4 we showed that quarks are confined vertically, in the 2-direction, 
but not horizontally, in the 1-direction, at the zeroth order of the weak- 
coupling expansion. To see what happens beyond this order, it is necessary 
to examine the quark-anti-quark potential in perturbation theory. This is 
very straightforward to do. 

If a quark is located at u 1 , u 2 , and an anti-quark is located at v 1 , u 2 with 
v 1 > u 1 , the electric- field operator in the 1-direction is given by 

l^x 1 ^ 2 ) = q{u l , u 2 )8 x i> u i8 x 2 u 2 - q(v l ,u 2 )8 x i> v i8 x 2 u i 

X 1 

- Y,(V 2 -l 2 )(y\x 2 ). (6.1) 
y l =o 

The constraint ()3.3j) is unmodified. Thus, the unperturbed states and ener- 
gies are unaffected by these two charges, as we claimed in Section 4. However, 
to first order in perturbation theory, there is a new contribution to Eq equal 
to 

AE^ = kCnIv 1 - u 1 ] , (6.2) 

where Cn is the smallest eigenvalue of the Casimir of SU(iV), q 2 = Cat 11, 
by ([2.5)1 . Thus, to first order in perturbation theory, the horizontal string 
tension is ^Cjv- What is especially remarkable about this result is that we 
can see clearly an electric string forming along the shortest path connecting 
the two quarks. 

What is happening physically is that the vacuum remains undisturbed 
by the charges and prevents the penetration of electric flux. To this low 
order of perturbation theory, we have a cost of at least mi to excite the 
chiral model at x 2 and x 2 — a. Thus there is a string tension equal to the 
(1 + 1) -dimensional string tension through an electric Meissner effect. We 
do not have to appeal to the condensation of some kind of magnetic charge 



13 



to make this interpretation. At higher orders of perturbation theory, the 
string of electric flux can presumably fluctuate; these corrections are needed 

„2 

to reliably set k = 

7 Conclusions 

In this paper, we have shown that lattice gauge theories in two space and 
one time dimension confine charges, through an anisotropic weak-coupling 
expansion. Though we cannot exactly evaluate the terms in this expansion, 
by just using some general knowledge of the chiral nonlinear sigma models, 
we can make precise statements about these terms. 

The astute reader may wonder if the methods developed here can work 
for the oldest known example of non-trivial confinement: lattice compact 
(QED)2+i PQ. The answer is that they do not. In this Abelian gauge theory, 
we would expand about the states of the U(l) nonlinear sigma model. This 
model has a massless phase at weak coupling, so we would not obtain vertical 
confinement and area-law behavior of the space-like Wilson loop. In fact, 
our perturbation method makes no sense at all for lattice (QED) 2 +i. The 
reason is that correlations of the operator l{x x , x 2 ) — l{x l , x 2 — a) (the adjoint- 
representation covariant derivative is simply the ordinary lattice derivative) 
do not fall off sufficiently fast to make ^2 ^1 directly proportional to the 
volume. The infrared divergence in the vacuum energy, which concerned 
us so much, really happens in the Abelian theory. This divergence is not 
real, but is an artifact of our methods. Our weak-coupling expansion seems 
peculiarly suited to non-Abelian theories in this regard. 

We have assumed that a mass gap exists in the (l + l)-dimensional SU(N) 
chiral model. At strong coupling, this can be proved rigorously with a cluster 
expansion, in the Euclidean lattice formulation. Perhaps a fully rigorous 
proof can be made of confinement with go large, but k small. 

Of the questions raised by our analysis we think that six stand out as 
important. We suspect, however, that only the first, second and possibly the 
third can be answered in the near future. 

The first and probably easiest important question is whether the infrared 
finiteness of our perturbation series exists beyond the second order. We hope 
to be able to settle this issue soon. If settled affirmatively, the series probably 
does not converge, but may be Borel summable. 

The second question is whether adjoint matter is confined for finite N. 
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This is certainly happening at first order in the horizontal direction. We 
believe that this property will disappear at higher orders. 

The third question is raised by the fact that our mass scales are set by 
()4. ip . with one exception (the horizontal string tension). All these quantities 
are non-zero for any positive value of a. We believe, however, that we should 
still have a mass gap and gap confinement as a — > 0, provided the continuum 
coupling constant go/(y/a) is kept fixed. Our vertical string tension, found in 
Section 4 is too small, and our horizontal string tension found in Section 6 is 
known only for small k. These numbers should both be proportional to g^/a, 
the square of the continuum coupling constant. Perhaps this difficulty can be 
removed by resummation of the perturbation series or by a renormalization- 
group argument. 

The fourth question is whether we can do a better job of calculating 
energies and states. Perhaps we could accomplish this, if Bethe's Ansatz for 
the chiral model could be carried out in a formalism where both the left- 
and right-handed SU(iV) symmetries are manifest in the Hamiltonian. In 
the work of Polyakov and Wiegmann [H] only one of these is manifest; the 
other appears in the S*-matrix, but its interpretation is obscure. If a version 
of Bethe's Ansatz with both symmetries manifest can be found, there is the 
possibility of a better understanding of the (2 + l)-dimensional gauge theory. 
One could use whatever regularization is most expedient for diagonalising the 
Hamiltonian, instead of the lattice. It may be a long time before this question 
can be seriously addressed (perhaps never). We believe a more likely path 
to success is to expand some version of the axial-gauge Hamiltonian about a 
system of (l + l)-dimensional field theories other than chiral sigma models. It 
would be a stroke of good luck, to have an expansion about exactly solvable 
field theories where the symmetries are easy to understand. 

The fifth question is whether our results can be understood in the context 
of condensation of magnetic charge. If a picture of condensing composite 
operators could work in the (1 + l)-dimensional chiral models (no one has 
succeeded in showing this), then operators defined on sets of points of one 
dimension higher should be important for confinement in (2 + 1) dimensions. 

The last and most important question is whether (QCD) 3+1 could be 
studied by our methods. This is, we hope to no one's surprise, a much 
harder problem. A lattice gauge theory in (2 + 1) dimensions is particularly 
amenable to the methods discussed here, because if the square of electric 
field in the 1-direction is dropped from the Hamiltonian, it easily breaks 
apart into (1 + l)-dimensional Hamiltonians we know a lot about. This 



15 



does not happen in 3 + 1 dimensions. The Hamiltonian breaks into (2 + 1)- 
dimensional Hamiltonians with both gauge fields and matter in the adjoint 
representation. These models are probably not even renormalizable, but 
seem worthy of investigation. 

Acknowledgements 

I thank Mike Creutz, John Kogut, Herbert Neuberger, and Tilo Wettig for 
organizing the workshop "Modern Challenges for Lattice Field Theory" , at 
the Kavli Institute for Theoretical Physics, where this work was completed 
amid the discussion of many interesting topics. I would also like to thank 
David Adams for discussions and a careful reading of the manuscript. 

This research was supported in part by the National Science Foundation 
under Grant No. PHY99-07949. It was also supported in part by a grant 
from the PSC-CUNY. 

References 

[1] A.M. Polyakov, Phys. Lett. B59 (1975) 82. 

[2] A.M. Polyakov, Nucl. Phys. B120 (1977) 429. 

[3] R.P Feynman, Nucl. Phys. B188 (1981) 479. 

[4] P. Orland and G.W. Semenoff, Nucl. Phys. B576 (2000) 627. 

[5] D. Karabali and V.P Nair, Nucl. Phys. B464 (1996) 135; Phys. Lett. 
B379 (1996) 141; D. Karabali, C. Kim and V.P. Nair B524 (1998) 661; 
Phys. Lett. B434 (1998) 103; Nucl. Phys. B566 (2000) 331, Phys. Rev. 
D64 (2001) 025011. 

[6] P. Orland, Phys. Rev. D70 (2004) 045014. 

[7] A.M. Polyakov, Phys. Lett. B59 (1975) 79; E. Brezin, J. Zinn- Justin and 
J.C. Le Guillou, Phys, Rev. D14 (1976) 2615; E. Brezin, J. Zinn-Justin, 
Phys Rev. B14 (1976) 3110; W.A. Bardeen, B.W. Lee and R.E. Schrock, 
Phys. Rev. D14 (1976) 985; Y.Y. Goldschmidt and E. Witten, Phys. 
Lett. 91B (1980) 392; A.B. Zamolodchikov and ALB. Zamolodchikov, 



16 



Ann. Phys. 120 (1979) 253; P. Wiegmann, Phys. Lett. B152 (1985) 
209. 

[8] A.M. Polyakov and P. Wiegmann, Phys. Lett. B131 (1983) 121; P. 
Wiegmann, Phys. Lett. B141 (1984) 217; L.D. Fadeev and N.Yu. 
Reshetikhin, Annal. Phys. 167 (1986) 227. 

[9] E. Abdalla, M.C.B. Abdalla and A. Lima-Santos, Phys. Lett. B140 
(1984) 71; P. Wiegmann; Phys. Lett. B142 (1984) 173. 

[10] S. Mandelstam, Bull. Am. Phys. Soc. 22 (1977) 541; Phys. Rev. D19 
(1979) 2391. 

[11] J. Kogut, Rev. Mod. Phys. 51 (1979) 659; M. Creutz, Quarks, Gluons 
and Lattices, Cambridge University Press (1983). 

[12] M.B. Halpern, Phys. Rev. D19 (1979) 517. 

[13] G.G. Batrouni, Nucl. Phys. B208 (1982) 12; M.B. Halpern and G.G. 
Batrouni, Phys. Rev. D30 (1984) 1782. 

[14] S.G. Rajeev, Phys. Lett. B212 (1988) 203; E. Langmann and GW. 
Semenoff; Phys. Lett. B296 (1992) 117; J. Hetrick, Int. J. Mod. Phys., 
A9 (1994) 3153; V. Miiller and W. Riihl, Nucl. Phys. B230 (1984) 
49; J.B. Bronzan, Phys. Rev. D31 (1985) 2020; J.B. Bronzan and T.E. 
Vaughan, Phys. Rev. D47 (1993) 3543. 

[15] A. McKane and M. Stone, Nucl. Phys. B163 (1980) 169. 



17 



